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Abstract 

We deal with a planar random flight {(X(t),Y(t)),Q < t < T} 
observed at n + 1 equidistant times U — iA n ,i — 0, l,...,n. The 
aim of this paper is to estimate the unknown value of the parameter 
A, the underlying rate of the Poisson process. The planar random 
flights are not markovian, then we use an alternative argument to 
derive a pseudo-maximum likelihood estimator A of the parameter A. 
We consider two different types of asymptotic schemes and show the 
consistency, the asymptotic normality and efficiency of the estimator 
proposed. A Monte Carlo analysis for small sample size n permits us 
to analyze the empirical performance of A. 

A different approach permits us to introduce an alternative estima- 
tor of A which is consistent, asymptotically normal and asymptotically 
efficient without the request of other assumptions. 

Keywords: asymptotic efficiency, discretely observed process, planar 
random flight, inference for stochastic process. 

1 Introduction 

Diffusion processes play a central role in the theory of stochastic processes. 
However these models do not give a realistic description of the real move- 
ments because the velocity is infinite and their sample paths are nowhere 
differentiable. For this reason in literature have been proposed alternative 
processes with finite velocity. The first model of this type, introduced by 
Goldstein (1951) and Kac (1974), is the telegraph process which describes 
the motion of a particle on the real line. 

The planar random flights are a possible extension in R 2 of the telegraph 
process. We consider the motion in the plane of a particle starting at arbi- 
trary point (xQ,yo), moving with constant velocity c and taking directions 
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uniformly distributed in (0, 2ir]. The changes of direction are governed by a 
homogenous Poisson process with parameter A > 0. Let N(t) be the number 
of Poisson events in the interval [0, t], the position at time t > of a particle 
performing a random flight is 

N(t)+1 

X{t)=XQ + c ( s j ~ s j-i) cos Qji 

3=1 
N(t)+1 

Y{t) = y + c ^2 ( s 3 ~ Sj-i)sin6»j, (1.1) 
j=i 

where 9j are independent random variables uniformly distributed in (0, 2ir], 
while Sj,j = 1, ...,n are the instants at which Poisson events occur (so = 
and SjV(f)+l = *)■ 

The distribution p(x, y; t) of (X(t),Y(t)) is concentrated on the disc 

S 2 ct = {(x,y) : (x - x ) 2 + (y - yof < c 2 t 2 }. 

If the initial direction is maintained until time t, the probability den- 
sity p(x,y;t) possesses a singular component, otherwise the distribution lies 
inside S 2 t . 

Random flights in R 2 have been studied in Stadje (1987), Masoliver et al. 
(1993), Kolesnik and Orsingher (2005). De Gregorio and Orsingher (2007) 
analyze random flights in M. d , d > 2, and derive their explicit distribution in 
the four-dimensional space. 

The only references about the statistical inference of random motion at 
finite velocity consider estimation problem for the telegraph process. Yao 
(1985) estimates the state of the telegraph process under white noise per- 
turbation and studies performance of nonlinear filters. Iacus (2001) consid- 
ers the estimation of the parameter 6 of a non-constant rate Xe(t). More 
recently, De Gregorio and Iacus (2006) introduce a pseudo-maximum like- 
lihood estimator and a least squares estimator for the parameter A when 
the sample paths of the telegraph process are observed only at equidistant 
discrete times. The authors also analyze the same statistical problem for 
a geometric telegraph process particularly interesting in view of financial 
applications. For a telegraph process observed at discrete times, Iacus and 
Yoshida (2006) study the asymptotic (i.e. the mesh decreases to zero and the 
horizon interval tends to infinity) properties of two moment estimators and 
propose also an estimator consistent, asymptotically normal and efficient. 

The aim of this paper is the estimation of the parameter A when the 
process {(X(t),Y(t)),0 < t < T} is observed at n + 1 equidistant times 
= to < t\ < .. < t n = T, where U = iA n = iA, i = 0, 1, n. We consider 
two types of asymptotic framework: 

1) A n — > and nA n = T — > oo as n — > oo; 
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2) n — ► oo with A n fixed. 

The statistical problem is interesting because the planar random flights 
seem to be useful for ecology and biology applications. In fact, Holmes 
(1993) and Holmes et al. (1994) consider these models to represent the 
displacements of the animals and microorganisms on a surface. 

We note that when the planar random flight is observed continuously, 
then N(T)/T is the optimal estimator of the parameter A and our statistical 
problem is equivalent to the estimate of a whole Poisson process on [0, T] 
(see Kutoyants (1998)). 

The process {(X(t),Y(t)),0 < t < T} is not markovian. Hence, it is not 
possible to explicit the likelihood function of the points observed as product 
of the transition densities. This fact implies that we can not use the tools 
developed for the diffusion processes (see Sorensen (1997) and Soresen (2004) 
for an account of these estimation methods). 

The main idea of this paper is to consider the points 

(X(iA n ) - X((i - l)A n ), Y(iA n ) - Y((i - l)A n )) 

as n independent copies of a random flight up to time A n (which is untrue) . 
In this manner we can build an estimating function from which it is possible 
to derive a pseudo-maximum likelihood estimator. 

The paper is organized as follows. In Section 2 we describe the random 
motion considered here and present some results concerning the process 
y/X 2 (t) + Y 2 (t) (for example the moments). In Section 3 we introduce a 
pseudo-likelihood function L n {\) and propose the following estimator 

A n = argmaxL n (A). (1.2) 

Under the asymptotic scheme 1) the estimator A n is asymptotically normal 
and efficient. Alternatively under the same asymptotic scheme , we present 
an estimator asymptotically efficient without additional hypotheses. By 
considering the second asymptotic framework in Section 4, we can study the 
convergence of the estimator A n by means of the pseudo-likelihood ratio. In 
the last section, we analyze the empirical performance of A n by means of a 
Monte Carlo analysis. 

2 Planar random flights: description and some re- 
sults 

We consider a particle starting at the arbitrary point (xo>£/o) of the plane 
M 2 , moving with constant finite speed c. The initial direction is a random 
variable 6 uniformly distributed in (0, 2tt]. The changes of direction are gov- 
erned by a homogeneous Poisson process with parameter A > 0. Therefore, 
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when a Poisson time occurs the particle takes a new direction uniformly 
distributed in (0, 2ir], independently from the previous one. 

We indicate the position of the particle at time t > with the stochastic 
process (X(t),Y(t)), which is called random flight in the plane. At time t 
the particle is located in the disc 

S 2 ct = {(x, y):[x- x ) 2 + (y - y ) 2 < c 2 t 2 }, (2.1) 



with probability 1. If no Poisson event occurs the particle reaches the circle 
dStf. = {(x, y) : (x — xo) 2 + (y — yo) 2 = c 2 t 2 }, with probability 

P{(X(t\Y(t))edS 2 ct } = e- xt . 



The remaining part of the distribution lies in the interior of (12. ip and 
represents the absolute continuous component of the distribution 

P {X(t) G dx, Y(t) G dy} . (2.2) 

We note that the random flights have trajectories which assume the form 
of broken lines where the single steps have random length and are uniformly 
oriented in (0, 2ir]. However, the total length for any sample paths at time 
t is ct. 

The density law of (X(t),Y(t)) (see Kolesnik and Orsingher (2005)) is 
equal to 



X e -At+fV c2 * 2 -(z-zo) 2 -(j/-J/o) 2 
yJcH z - (x - x ) z - {y - yo ) 

+ e -^5{c 2 t 2 -{x- xo) 2 - (y - y ) 2 ), (2.3) 
2nc 

with (x,y) G S 2 t and <5(-),l(-) representing respectively the Dirac's delta 
function and the indicator function. 

Now, we present some results concerning the following process 



R(t) = y/X*(t)+Y*(t), (2.4) 

i.e. the euclidean distance from the origin of the space IR 2 of the position 
reaches by the moving particle at time t > 0. 

The following theorem contains our first result. 

Theorem 2.1 The absolute continuous component of the process R(t),t > 
0, when R(0) = \/ Xq + y 2 , is equal to 



\ p2n ^y/c 2 t 2 -r' 2 -x 2 -y 2 +2r(x cos 0+2/0 sin 0) 

f R (r,t) = ^re- xt / =d9, 
27rc Jo v ° 2 * 2 ~~ r2 ~~ x o ~ Ua + 2r(xo cos 6 + yo sin 0) 

(2-5) 

mf/i <r < ct. Moreover, under the Kac condition (i.e. c, A —* oo in sitc/i 

2 

a way i/iai y — >• 1^, we /lave £/ia£ (|2.5p tends to the law of a standard Bessel 
process. 
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Proof. We start the proof observing that 

P {R(t) <r\R(0) = y/xl + y*} (2.6) 



e -\t+±y/cH2-(x-x ) 2 -(y-vo) 2 

zdxdy 



JJ{( Xt yy. x 2 +y 2< r 2 } ^/ c 2 t 2 - (x - x ) 2 -(y- yo) 2 

\ rr j'2-K g— At+-^ \J c 2 t 2 — {x$—p cos 9)' 2 — (yo — pcos d)" 2 

/ pdp / — dO 

27rc Jo Jo ^/c 2 t 2 - (x - p cos 6) 2 - (yo - p cos 6) 2 



^ f-r f-2-K e -At+|- v / 'c 2 t 2 - p 2 -x 2 -y 2 +2p(x cos6»+J/o sintf) 

- — / pdp / — =d9. 

27rc Jo Jo \J(?t 2 - p 2 - Xq - yQ + 2p(xo cos 9 + y sin 9) 



By differentiating the probability (|2.6p with respect to r, the density 
2.5p emerges. 

In order to prove the second part of the theorem, we rewrite the density 
2.5p in the following form 



r 2 +xi + y 2 -2r(xQ cos 9 + j/q sin 6) 

2vr „- Xt+Xt \l 1 ^Tp 



A f Z7T e 

fdr,t) = ^-r , =— = = = d6 



27TC J L r 2 +x 2 +y 2 - 2r(x cos 0+y sin 6) 
Ct\l 1 -375 



c 2 t 2 

r 2 +x 2 ) + y 2 — 2r(xQ COS0+&Q sing) 



A ^ e 2 ?< 



r 1 ....... . ^9. (2.7) 



2lTC 2 t Jo /~. r 2 +x 2 +y 2 — 2r(xo cos s i n 
V 1 



In the last step we have used the expansion y/l — w = 1 — — ^ — which 
is absolutely convergent for |io| < 1. 

From ()2.7p . under the Kac condition, we obtain the following limit 



2?r 2 ^2_,_„2 



r / _r£ *(]+«() rt^QCOsfl + j/Qsinfl) 

hm f R (r,t) = — - / e 2t e 2t e * ^ 

A,c^oo lilt 
c 2 /A^l 







= -e 2t e » I — S£ — ^ , (2.8 

by means of the well-known integral representation of the Bessel function 

f2w 



I (x^a 2 + f3 2 ) = — I n e^ acos9+l3sin 
27r Jo 



Expression (12. 8|) represents the density function of Bessel process B 2 (t) + B^t), 
where B\ and B2 are two independent standard Brownian motion. □ 
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Prom Theorem 2.1 we note that for (xq,uq) = (0,0), the complete dis- 
tribution of R(t) becomes 

A rexp{-At+|Vc 2 t 2 -r 2 } 
P*(r,t) = ^2372 ^<r<^ ( 2 - 9 ) 

+ r -^5(c 2 t 2 -r 2 ). 
ct 

We note that the density (|2.9p coincides with formula (7) in Kolesnik 
and Orsingher (2005), when we ignore the angular component. 

By taking into account the probability law (12. 9p . we are able to derive 
the moments of R(t). 

Theorem 2.2 Let (xo,yo) = (0,0) and p > 1, we have that 

E^(t) = (ctf e - A '|v^(A) 2 r(^±i-)/ 2 ± i (At) + l|. (2.10) 
Proof. In view of (|2.9p . we can write 

\ ret -\/c 2 t 2 -r 2 

VRP( t ) = -e- Xt / r? +1 J dr + {ctfe~ Xt . (2.11) 

c Jo Vc 2 t 2 - r 2 

Now, we work out the integral in (I2.11|) . Hence 

ct e ±Vc 2 t 2 -r 2 
r v +l — r dr 
V c 2 t 2 — r 2 

- £ I {^) k C rP+1{cH2 ~ r2)hfldr = {r = ct ^ ] 

1 /A\ fc (ctr+ fc+1 r(| + i)r(Y) 

^r(g + i)vl^V m =(* = *») 

b )^k\\2j r(Mi + | + i)r(|) 1 J 

V2 ) y 1 ^{2) 2m r( m) r (m + £±i + 1 

p+1 



v-r(^ + i)( c ^(-j /^(At). (2.12) 
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By inserting (|2~T2l) into ([2~PT]) we obtain the result (f27TOl) . 



□ 



Remark 2.1 We observe that 



A— >oo 



lim ER p (t) = 0. 



(2.13) 



In other words, if A grows to infinity the changes of direction increase 
and consequently the distance from the origin decreases. 

Remark 2.2 We derive from (|2.1U|) the mean value of R(t) 



In the particular case p = 2, we can write the square mean in terms of 
simple function. In fact, by means of the following relationship 



3 Parametric estimation for planar random nights 

We assume that the planar random flight {(X(t), Y(t)), < t < T}, with 
(X(0),Y(0)) = (0,0), is observed only at n + 1 equidistant discrete times 
= Iq < t\ < ... < t n = T, where ti = iA n = iA, i = 0,1,..., n. 
We use the following notation to simplify the formulas: (X(ti), Y(ti)) = 



The interest is the estimation of the parameter A whilst the velocity c 
is assumed to be known. In other words, we want to estimate the rate of 
change of a microorganism which performs a planar random flight, when we 
are able to observe its position only at discrete times. 

The estimation of c is an uninteresting problem. In fact, if in the interval 
(iA n _i,iA n ] there are not changes of direction, then (Xi — Xj_i) 2 + (Y{ — 
Yi-i) 2 = <? A\. If A n is suitable small, there is high probability of observing 
N(ti) — N(ti-i) = and c can be calculated without error. 

Analogously to the telegraph process, the random flights are not marko- 
vian. For this reason we cannot write the explicit likelihood of the process in 
the form of product of transition densities as well as for diffusion processes. 
Therefore, we need an alternative argument in the spirit of the paper by De 
Gregorio and Iacus (2006). 



ER(t) = cte- xt {^h{\t) + 1}. 



(2.14) 




we get that 




(X{iA n ),Y(iA n )) = (X u Yi). 
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n 

-10 



10 



X(t) 



Figure 1: Discrete time sampling of the planar random flight. For this 
sample path n = 7 and n + = 6. 



We define a pseudo-likelihood function as follows. By taking into account 
the distribution (j2.3l) . we introduce the following data dependent function 



L n (X) 



L n (X\(X ,Y ),(X 1 ,Y 1 ),...,(X n ,Y n )) 

n 

11 p{ (Xt , Yi ) , A n ; (Xi_i , ) , U- 1) 
A exp{-AA n + ^y/u^} 



(3.1) 



i=l 
n 

n 

i=l 



-AA„ 



2vrc 



/Mr, 



L {«n,i>0} 



+ 



2vrc 



= 0) 



where u n ,i = u n ({X h Y { ), (X^Y^)) = c 2 A 2 n - (X { - X^) 2 - (Y, - Y^) 2 . 

The transition densities p((Xi,Y{), A n ; appearing in 

()3.ip represent the distribution of a random flight in M 2 , initially located 
at the point (JQ_i,Yi_i) at time which reaches the position (Xj,Yj) 
at the instant £j. The function (|3.ip is indeed the joint law of the points 
{Xi — Xi—i,Yi — Y-i)> which are considered as if they were n independent 
copies of the process (X(A n ), Y(A n )) (i.e. the process (X(t),Y(t)) up to 
time A n ). 

The pseudo-likelihood function (|3.ip is equivalent to 



MA) 



n 



A exp{-AA„ + ^Ju^j} 



2irc I AJ- I 27rc 

-AnA n A n+ exp I £ X^r=i 



(2vrc) r 



nr= + i 



(3.2) 



/Mr, 



where n + is the number of the planar random flights with at least one change 
of direction. 
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Remark 3.1 In the expression (13.20 . the factor f e 2 " j concerns 

the singular part of the densities p((Xi, Yi), A n ; (Xi—i, li_i), tj-i), while the 
product represents the absolutely continuous components of the distribu- 
tions of the random nights. Note that for increasing values of A, the abso- 
lutely continuous component of (13. 2p has a bigger weight than the singular 
component; viceversa for small values of A. Figure [T] shows how the two 
components of the function L n (X) emerge for this scheme of observation. 

To define our estimator, we borrow the approach based on estimating 
functions. Formula (|3.2p yields 

d 1 n+ n + 

F nW = -rrlogL n (\) =-uA n + -V^+^, (3.3) 
a\ c A 

i=i 

and the estimator is obtained by solving F n (A) = 0. 

By taking into account the function (|3.3|) we derive the following pseudo- 
maximum likelihood estimator for the parameter A 

A„ = argmaxL n (A) 



cn 



cnA n - Yd=l sl^r, 



(3.4) 



It is easy to see that ^-logL n (A) < and the uniqueness of the esti- 
mator (|3.4p holds. However A n is not a true maximum likelihood estimator 
because L n (X) is not a true likelihood function. In the section 4 we will 
analyze the empirical performance of the estimator (|3.4|) for small sample 
size. 

By considering the following asymptotic framework A n — > and nA n = 
T —* oo as n — > do, we provide the next result for the estimator (13. 4D . 

Theorem 3.1 Let A n — > and nA n = T — > cxd as n — > oo, then X n is 
consistent, asymptotically normal and efficient. 

Proof. If A n — > and nA n = T as n — > oo, we have that u n i — > and 
n + — > N(T), where N(T) represents the number of changes of direction 
occurred in the interval [0, T] when the whole trajectory is observed. Thus, 
A n tends to the maximum likelihood estimator of a homogeneous Poisson 
process 

^1. (3.5, 

Therefore, under the condition nA n = T — > cxd, the pseudo- maximum 
likelihood estimator is consistent, asymptotically normal and efficient (see 
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Kutoyants (1998)). □ 

We introduce an alternative estimator for the parameter A by means of 
the distances 

rn = y/fXi - X^f + {Yi - Yi-{f. 

By setting 



-y n 1 n 

Gn = ^ E Hvi^An} = ^ E 1 
i=l i=l 

we define the following unbiased estimator 



{N([u-i,ti))>l}' 



A„ = — log (1 - A n G n ) . (3.6) 

The advantage of A n is that we are able to derive the asymptotic prop- 
erties without assumptions on the points (Xj — Xj_i, 5^ — l^-i). 

Theorem 3.2 For nA„ = T — ► oo and A n — > as n — > oo the estimator 
(|3.6p is consistent, asymptotically normal and efficient 



V^A^(A n -A) 4jV(0,A). (3.7) 

Proof. We replace the steps contained in the proof presented by Iacus and 
Yoshida (2006) for the telegraph process. 

First of all we prove consistency and asymptotic normality of G n . We 
observe that 

1 — p—^An 1 

E(G») = An = A + -A 2 A n + o(Al) - A 

and consistency immediately follows. Now we show the asymptotic normal- 
ity of G n and consider to this scope the quantity 



U n = ^fnA n {G n -V(G n )) 
1 n 

= E [ 1 {%<ca„> - E(i {j?i<cAn })] 

1 n 

= ^Eiw.^n-d-e-^) 

n 

= £ 



a; 
=1 



where 



a; 



' l {J v([^ lit0 )>i } -(l-e- AA ") 



VnA n 
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It's clear that E(aij) = and E(U n ) = 0. Moreover 
Var(ai) = -^-Var(l {iV ([t i _ 1 ,t,))>i}) 



{E(l{7V([ il _ 1 , tl ))>i}) - (E(l {iV ([ ii _ lA) )>i})) 2 } 



nA r 
1 

1 

nA n 

~{A + (1)} 

n 



e -AA n _ g -2AA„ 



therefore 



i.e. 



Var(C/ n ) = A + o(l). (3.8) 
The variables «j are independent and the Lindeberg condition is true, 

E E {V([^iA))>i}«?}^ ' (3-9) 

i=l 

because for large n it holds true that laA < , . . From condition (13.91) 
follows that 

l7 n -i]V(0,A). (3.10) 
Finally, we can prove the asymptotic normality of X n . Since 

/H = --^log(l-™A n ), f'{w) ' 



n 1 - u;A n 

and 

A„ = /(G n ), A = /(E(G n )), 
then, by so-called <5-method, we have that 



V^A^-A) = ynK n {f{G n ) - f(E(G n ))) 



y/vAn{G n - E(G„))/'(A) + o p ( ^/nA n \G n - E(G n )\) 

1 

1 -~AA 



y/nA^(G n - E(G n ))- — + o p (v^\G n - E(G n )\), 



hence for nA n = T — » oo and A n — * as n — * oo, we obtain that 

\/^(A„-A)4 iV(0,A). 

□ 
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4 Large sample properties for the pseudo-maximum 
likelihood estimator 



To analyze the properties of the estimator (13, 4p as n — > oo with A n fixed 
(large sample scheme) , we use the tools of asymptotic theory of statistical 
estimation presented in Ibragimov and Has'minskii (1981). 

Let us assume the velocity c known and A G (Ai,A2) = with < 
Ai < A2 < 00. We need to introduce another hypothesis: the distance A n 
between two consecutive instants t{,i = 0, 1, n, is such that the following 
condition holds 

P A {(X(A n ),Y(A„)) e mtS 2 cA J = 1, (4.1) 

where int S 2 A = {(x, y) : x 2 + y 2 < c 2 A 2 }. 

In other words between the points (Xi_i,Yj_i) and (Xi,Yi),i = l,..,n, 
the planar random flights have at least one change of direction (or equiva- 
lently P\{N(iA n ) = 0} = 0, i = l,...,n). In general it is obvious that for 
increasing values of A, the minimum value of A n satisfying the condition 
(|4.ip decreases. 

Immediately, from (|4,ip follows that the singular part of (|3.2p vanishes. 
In fact we have that 



Ln(X) = n 



1=1 



A exp{-AA n + l^/n^l} 
2ttc i/u n .i 



t A \ n exp { -AnA n + | Y^=i ^i] 



(4.2) 



while the pseudo-maximum likelihood estimator (|3.4p becomes 



An = T =■ (4-3) 

We start our analysis observing that the Radon-Nikodym theorem yields 

A 



dP\ ( A \ exp (-At + -,/WTi) 
d/i \ 27rc / 

J4.4) 

where /z is the Lebesgue measure in the plane. Thus, we can indicate L n (X) 
as follows 

MA) = ^, («) 

where represents the joint probability distribution of n independent 
copies of a planar random flight up to time A n . 

It's appropriate to remark that we are presenting in this section results 
valid only for the parametric model {P^,A € 0}, i.e. the model deriving 
from the assumption of i.i.d. observations. 
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To simplify the formulas we write p((Xi, YJ), A n ; l^-i), U-\) = 

p{\). Our first result is the following theorem. 

Theorem 4.1 Let £ n be the experiment generated by n independent obser- 
vations of (X(A n ),Y(A n )). Then £ n is regular with Fisher's information 
equal to 

n 

InW = (4-6) 

Proof. By considering the definition of regular experiment presented in 
Ibragimov and Has'minskii (1981), page 65, we must prove that \/p(X) is 
differentiable in L2 (the space of the square integrable functions) with con- 
tinuous derivative in L2 



m = (-An + - c v^+\ ) < 1-7) 



By setting g(X) = y/p(X) we get that 

If (g{\ + h)-g(X)-H(X)fdxdy 



(4.8) 



Now, by observing that 

-^+^ v ^7+io g yA±F _ h ( A 1 — r 



l + ^(-A n + ^^H—+^j +o(h), (4.9) 



we obtain 



I ! (g(X + h)- g(X) - h^{X)fdxdy = o(\h\ 2 ). 

The continuity of ip(X) is shown by means of the dominated convergence 
theorem. 

To complete the proof we verify that £ n possesses finite Fisher's infor- 
mation I n (A) for any A € G. Clearly I n (X) = nl(X), where I(X) represents 
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Fisher's information of a single experiment. Thus, we can write 



/(A) 



7/ 2 |V;(A)|: 



dxdy 



A 
2vrc 

A 
2vrc 



-AA n + -./u^T 



Tl.t 



- - A n + -^/u^ dxdy 
A c v ' 



-AA„ 



r._> 



III 



n.i 



1 . \ * U ni 2/1 



cA r , 



— / « 



c V A 



2tt pe -AA„+A v ^Af^2 



\A 2 A 2 -p 2 

*n P 



2 cPAl-o 2 2/1 



(4.10) 



+ c A 



where in the last step we have used the transformation in polar coordinates 
x = xq + p cos 8,y = yo + p sin 6. 

To obtain the explicit value of (|4.10p we calculate the following three 
integrals 



1\ 



A /l 



2vrc V A 



A /l 



A, 



2 r cA 



-It -A, 



c V A 



5-* 

5-* 



n 


[■2tt 











rcA n 








Jo 





V^A 2 -/; 2 



-AA, / A^/ c 2 A 2_ p2 



1 - e 



p=cA„ 



p=0 



-AA r 
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J 3 = J^dpfj dOpe-^^^ 



2A ( l A 

c 2 " ( A _ 



,-AA„ 



pe 



^ A l-P 2 dp =(z= \/c 2 A 2 - p2) 



cA r , 



-AA n 



2 A 2 



1 



2I--A, 



2=cA„ 2 

2=0 C 

2 A. 



A n e" AA " 2A n e AA « - ±e 



A 





2=cA n 



e^dz 



z=0 



A, 



A 



1 - e 



-AA„^ 



Putting together X\ , %i , T3 we have that 



1(A) 



1 - e 



-AA n 



h A 2 

A 2 n 



By assumption P^{A r (A n ) = 0} = e 



A ^ (l " e" AA " (1 + A 2 A 2 



-AA, 



the result P~5]) follows. □ 



Fisher's information plays a central role in the Cramer-Rao inequality 
and more in general in the parametric inference. Let E^(-) be the expecta- 
tion with respect to the probability measure P™. For any estimators of the 
parameter A, we have the next result. 
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Theorem 4.2 LetT n be an arbitrary estimator of X such that E™|T n | 2 < oo 
for any A > 0. Then 

6(A) = E n x T n - A (4.11) 

is differentiable respect to A in L2. Moreover, the following Cramer-Rao 
inequality holds 



(I 1 db(X) n2 

Bl(T n - A) 2 > ' ,/A 



/n(A) 



+ 6 2 (A). 



(4.12) 



Proof. We note that 



ElT n = B n Xo \T n [y) exp [ -(A - X )nA n + 



A-A 



i=l 



and show that E-?T is differentiable and the equality 



i E " T »= E M T »(i-" A 



1^, \\ 



(4.13) 

(4.14) 



holds in L2. 

For this purpose we interpret ^E^T n as the right-hand side of equation 
(14.141) . It is not difficult to see that 



En rp -prirp t rpnT' 



1-h 



n 



nA 

A c 



1 n 

>n ~t~ — ^ ^ \[^n. 



i=l 



< E^|T n | 2 E^ 



V 1=1 



+ o(fc), 



where in the last step we have used the Cauchy-Schwarz inequality. 
By inserting the equality 

e -hnA n + ± £» =1 v^T+nlogd+VA) = 1 + /j ^ _ n A n + ~ ]C 

into (|4.15p . we can conclude that 6(A) is differentiable in L2-sense. 

The validity of the inequality (I4.12D follows by standard arguments. □ 

Remark 3.3 By taking into account an unbiased estimator T n of the 
parameter A, from (|4.12p we get that 



E™(T n - A) 2 > -. 



(4.16) 
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It's well-known that the Cramer-Rao doesn't give a good definition of 
asymptotic efficiency, because the limit variance may not coincide with the 
variance of the limiting distribution. Therefore to investigate the asymptotic 
properties of the estimator A n as n — > oo and A n fixed, we reduce our 
problem to the study of the normalized pseudo-likelihood ratio 

dP n 

ZnAz) = ~^ L = (4-17) 

A 

exp iy^r~ v^M - f(n)zA n + log ^ J J 

( (p(n)z ^ / 

I — - — 2^ V^M _ <f>{n)nzA n + ralog ^ 



exp 



A + ip(n)z 



A 



where <p(n) = <p(n, A) = (J n (A)) 1 / 2 . The function (|4.17p takes values in the 
following set 

U n A = < z : A + - = eel. 

It's well-known that ^ n> A (deriving from an i.i.d. observation scheme) 
admits the representation 

with , 1 X^ILi -^(0) 1) an< i </>n(z, A) — > in probability as n — > 

oo; i.e. is locally asymptotically normal (LAN). 

For the function Z n \ we have the next useful Lemma. 

Lemma 4.1 Let K be a compact subset of Q. We have that: 

i) for some constant a = a(K), B = B{K) 



sup sup \z — v | 2 E™ 

XeK \z\<R,\v\<R 

with z, v G U n± \\ 



Zllliz) ~ Z'jliv) 2 < 5(1 + ^), (4.19) 



/or any z G f7 njA 



z) < e 



-c|*| a 



(4.20) 



where c > 0. 
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Proof, i) Following the proof of Lemma 1.1, section III in Ibragimov- 
Has'minskii (1981) we get that 

2 



< 



4/2, 



(4.21) 



Since 



for z > follows 



A 

A+^ 



1+ A 



-2 



I n (A + y(n)(2 + s(u - z))))ds 
*(* + *) f A V 

/(a) ^a+*; ' 

< 1. If z < we can see that 
:l + (-2)i+ofi) <3 + o(~ 



z — v 



A 



Therefore set B = 3 + o (f) the inequality 



sup sup 

Aee |z|<i?,,A+^ee 



in (A + z) 



4(A) 



< 5(1 + 



(4.22) 



holds. 

In view of the relationships (|4.22p and (|4.21j) the proof of the inequality 
(I4.19p is concluded. 

ii) The function d^Jp(X) / d\ is differentiable in L2, then 

2 



S' 2 



\y/p(X + h)-y/pJXj\ 2 dxdy 



h 2 



°cA r , 



d_ 



y/p{\) dxdy + o(\hf 



-I(\) + o(\h\ 2 ). 



By taking into account that 



< inf I(X) < sup /(A) < 00, 



we have immediately that 



\y/p(\ + h)- y/pJX)\ 2 dxdy > 0. 
From ()4.23p . we derive the inequality 



S cA„ 



(4.23) 



inf inf 

\£K {h.X+heS} J J S 2 



Wp{\ + h)- ^pJX)\ 2 dxdy > 



a\h\' 



l + |/i| 2 ' 



a > 0, 



and Lemma 5.3, Chapter I, in Ibragimov and Has'minskii (1981) permits us 
to obtain the condition (|4.2U|) . □ 



Finally, we are able to present the main result of this section. 
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Theorem 4.3 Let K be a compact subset of 0. Then the estimator X n , 
defined in (|4.3p . uniformly in A G K : 

• is consistent; 

• converges in distribution as follows 

V^(A)(A n - A) 4 iV(0, 1); (4.24) 

• /ias moments such that 



lim E^|/UA)(A„ - A)P = E|e| 7 , (4.25) 
where 7 > and £ ~ N(0, 1). 

Proof. In accordance with the Theorem 1.1 Chapter III, in Ibragimov and 
Has'minskii (1981), we prove that the four conditions are satisfied. 

The probability measure is uniformly local asymptotic normal, while 
it's easy to see that 

lim sup <p 2 (n, A) = 0. 

The validity of Lemma 14.11 concludes the proof. □ 

The Theorem 14.31 yields the Hajek-Le Cam asymptotic efficiency of the 
estimator A n with respect to a quadratic loss function. In fact, we have that 



lim lim sup v / /„(A)(A ri , - A) =1. (4.26) 

S^On^oo |A-A |<<5 

5 Monte Carlo analysis 

We analyze the empirical performance of the pseudo-maximum likelihood 
estimator A n by means of a Monte Carlo analysis with n < 00 fixed. We 
simulate 10000 sample paths of the planar random flights in the interval 
[0,T], with T = 500, for different values of A and c = 1. For any trajec- 
tories we have sampled n = 200, 300, 500, 1000 values subsequently used to 
estimate the unknown parameter A. 

The results have been reported in the Table [TJ Furthermore in the Table 
[T] there is a column y / MSE(A) derived as follows 



VMSE(A) 



\ 



1 N 

-]T(A n -A)2, (5.1) 



4 = 1 

where iV = 10000 is the number of simulations. 
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It emerges, as expected, that the mean square error tends to zero when 
the sample size increases. Furthermore, it is clear that the true value of the 
parameter A and the mean square error are correlated. In fact, for fixed 
n, as the more A increases the more Poisson events remain hidden to the 
observer. The bias assumes small values for all the cases considered and is 
constantly equal to 0.002 for A = 0.1, 0.25, 0.5, 0.75. 
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A 


Bias 




min A n 


max A n 


n 


0.10 


0.002 


0.015 


0.05 


0.15 


200 




n nnn 

0.002 


0.015 


n n cr 

0.05 


0.16 


onn 

300 




0.002 


0.015 


0.05 


0.15 


500 




0.002 


0.014 


0.06 


0.15 


1000 


/ \ or 

0.25 


0.002 


0.026 


0.17 


0.37 


200 




0.002 


r\ no V 

0.025 


n 1 7 

0.17 


0.35 


onn 

300 




0.002 


0.024 


0.17 


0.35 


500 




0.002 


0.023 


0.17 


0.34 


1000 


0.50 


0.001 


0.042 


0.37 


0.67 


200 




0.002 


0.038 


0.36 


0.65 


300 




n nnn 

0.002 


r\ no V 

0.035 


0.36 


0.65 


mn 

500 




0.002 


0.033 


0.36 


0.63 


1000 


0.75 


—0.000 


0.057 


0.56 


1.05 


200 




0.001 


0.051 


0.56 


0.98 


300 




0.002 


0.046 


0.60 


0.99 


500 




0.002 


0.042 


0.62 


0.93 


1000 


1.00 


—0.004 


0.073 


0.76 


1.28 


200 




—0.001 


0.064 


0.76 


1.29 


300 




0.001 


0.056 


0.81 


1.26 


500 




0.002 


0.050 


0.82 


1.18 


1000 


1.50 


-0.013 


0.106 


1.15 


1.98 


200 




-0.003 


0.090 


1.19 


1.92 


300 




0.001 


0.076 


1.22 


1.78 


500 




0.001 


0.066 


1.26 


1.78 


1000 


2.00 


-0.031 


0.141 


1.49 


2.53 


200 




-0.010 


0.117 


1.57 


2.61 


300 




0.000 


0.097 


1.67 


2.41 


500 




0.001 


0.080 


1.69 


2.29 


1000 



Table 1: Empirical performance of the estimator A n defined in (|4.3|) for 
different values of the parameter A and different sample sizes. The velocity 
c assumes value 1. The time horizon T is equal to 500. The results have 
been obtained on 10000 Monte Carlo sample paths of the planar random 
flights. 
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